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ABSTRACT 

A strong outburst of the October Draconid meteor shower was predicted for October 8, 
P^ . 2011. Here we present the observations obtained by the Canadian Meteor Orbit Radar 

["tI ' (CMOR) during the 2011 outburst. CMOR recorded 61 multi-station Draconid echoes 

and 179 single-station overdense Draconid echoes (covering the magnitude range of 
V6 < My s^ -f7) between 16-20h UT on October 8, 2011. The mean radiant for 
the outburst was determined to be ag = 261.9° ± 0.3°, 5g = -^55.3° ± 0.3° (J2000) 
O . from observations of the underdense multi-station echoes. This radiant location agrees 

pi ' with model predictions to ~ 1°. The determined geocentric velocity was found to be 

^2 ' ~ 10 — 15% lower than the model value (17.0 — 19.1 km • s~^ versus 20.4 km • s~^), a 

discrepancy we attribute to undercorrection for atmospheric deceleration of low density 
Draconid meteoroids as well as to poor radar radiant geometry during the outburst 
peak. The mass index at the time of the outburst was determined to be ~ 1.75 using 
the amplitude distribution of underdense echoes, in general agreement with the value 
of ~ 1.72 found using the diffusion-limited durations of overdense Draconid echoes. 
\l ', The relative flux derived from overdense echo counts showed a similar variation to the 

meteor rate derived from visual observations. We were unable to measure the peak 
flux due to the high elevation of the radiant (and hence low elevation of specular 
Tij" , Draconid echoes). Using the observed speed and electron line density measured by 

^D ' CMOR for all underdense Draconid echoes as a function of height as a constraint, we 

CO I have applied the ablation model developed by Campbell-Brown & Koschny (2004]) • 

From these model comparisons, we find that Draconid meteoroids at radar sizes are 
consistent with a fixed grain number Ugrain = 100 and a variable grain mass rugrain 
between 2 x 10~^ kg to 5 x 10^^ kg, with bulk and grain density of 300 kg • m~^ and 
3 000 kg • m~'^ , respectively. One particular Draconid underdense echo displayed well- 
^ I defined Fresnel amplitude oscillations at four stations. The internal synchronization 

- - allowing us to measure absolute length as a function of time by combining the absolute 

timing offsets between stations. This event showed clear deceleration and modeling 
suggests that the number of grains for this meteoroid was of the order of 1 000 with 
grain masses between 10^^° to 10~^ kg, and a total mass of 2 x 10^^ kg. 

Key words: meteors, meteoroids, comets: individual: 21P/Giacobini-Zinner. 



1 INTRODUCTION cal history. It is also among the very first (if not the first) 

meteor shower to be predicted on the basis of the parent 

1.1 A Brief History of Draconids , , ■ -j- j- r? iu: ui • ^ iu u u ■ 

■^ comet s proximity to Earth s orbit prior to the shower being 

Although quiet in most years, the October Draconid me- widely observed. 

teor shower (sometimes shorten to "Draconids", and also 

referred to in older literatures as the "Giacobinids" , after its 21P/Giacobini-Zinner was discovered by Michel Gia- 

parent body 21P/Giacobini-Zinner) has produced several of cobini in 1900, before any modern observations of the Dra- 

the most spectacular meteor storms in modern astronomi- conids. Due to limited observations, the exact orbit of the 

comet was un known until Ernst Zinner recovered it in 1913 

IjKronkl 120081 '). The orbital period was determined to be 

* E-mail: qye22@uwo.ca ~ 6.5 yr, identifying the comet as a Jupiter-family comet. 
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a major comet family that occasionally experience strong 
perturbations from Jupiter. 

The possibility of a mete or show e r gen erated by the 
comet w a s first proposed by iDavidsonI (JIOISJ), years before 
iDenningj (| 19261 ') first observed the meteor shower in either 
1920 or 1926 (Denning only saw 5 meteors in 1920, therefore 
the linkage is questionable). The maximu m zenith hourl y 
rate (ZHR) in 1926 was reported to be 17 (|Denningjll92^ . 
Before the historic Draconid meteor storm in 1933, no other 
definitive observations of the shower had been reported. 

The 1933 Draconid meteor storm was largely a sur- 
prise as no specific prediction suggested an unusually in- 
tense shower that year. The meteor rate reported by Eu- 
ropean and A merican observers was between 2 400 hr~ ^ 
and 5 400 hr'^ (jWatsonlll93ilNiiland fc van der Biltlll935h . 
while more recent estimates have suggested 10 000 ± 2 000 
IjJenn iskens 1995). The apparent radiant was measured to 
be near a^ = 266°, 5^ = -h59° bv iNiiland fc van der Bii^ 



1935 ). with a relatively large diameter (~ 10°). iFishen 



19341 ) suggested as many as six possible historical sight- 



ings before 1920, with the earliest sighting dating back to 
Chinese observations in 585 A.D ., although later studies do 
not support that hypothesis fe.g. llmoto k, Hasegawalll958l lH 

A Draconid outburst at 21P/ Giacobini - Zinne r's 1939- 
1940 apparition was predicted by IWatsonI l|l939r ). but no 
activity was observed in either year. However, Watson made 
a successful out burst predictio n for 21P/Giacobini-Zinner's 
1946 apparition (|Watsonlll946l ). The Draconid meteor storm 
again favored American and European observers, with max- 
imum activity at around l Oh UT on Octobe r 10, 1946 and 
a ZHR of 12 000 ± 3 000 jjenn iskcns 1995). It was also a 
landmark in meteor studies, not only because it was the first 
time that radars were deployed to observe a meteor sh ower 
l|LoveU et ai]|l947l : lHev et al.ll 19471 : IStewart et "aLlll947l '). but 
also because during this meteor storm many precise photo- 
grap hic records of brigh t Draconid meteors were obtained 
(e.g. Ijacchia et al.lll95Cll ). 

The 1946 Draconid meteor storm was the second and 
the last recorded Draconid meteor storm in the 20th cen- 
tury. 21P/Giacobini-Zinner's apparitions in 1985 and 1998 
were accompanied with sub-storm outbursts with a max- 
imum ZHR of ~ 700 iLindbladl [19871 : lArld Il998l ). Mod- 
erate activity was detected by radar in 1952, 1972, and 



^ The Chinese observations were quoted from Sui Shu (the Book 
of Sui), Vol. XXI, Sect. 3, which recorded that a meteor outburst 
took place on Sept. 23, 585 A.D., approximately corresponding 
to Aq = 202°in J2000 epoch. Most studies mentioning this event 
quote lBiotI 1 1184 ll ): "...hundreds of meteors fled in every direction." 
But there was one sentence from the original Sui Shu text omitted 
by Biot: "(astrologers) divined and said: 'small meteors fleeing 
in every direction means common people is moving around.' ", 
appar ently indicating the meteors were mostly faint. Ijenniskenj 
(120061 ) suggested a linkage to the Orionids; but given the fact that 
the moon was at its last quarter, and the Orionid radiant didn't 
rise until midnight, the event was unlikely to take place other 
than in the early evening as the moon would wash out most faint 
meteors. Since the observations were conducted in the capital of 
Sui, which was Chang'an (today's Xi'an); at a latitude of 34° 
N, the radiant of the Draconids would be high enough for the 
described viewing geometry ("fled in every direction"), therefore 
a linkage of the 585 A.D. observation to Draconids cannot be 
ruled out. 



2005 (iDavies fc Lovelll Il955l : iHughes fc ThompsonI Il973l : 
ICampbell-Brown et al.ll2006l ). but the associated visual ac- 
tivity was either unreported (195 2), very weak (1972), or 
moderate (2005) with ZHR-^ ^ 40 (JMalev fc Saulietid [l973 : 
ICampbell-Brown et al.l |2005| ). In other years, the hourly 
rates of Draconids were no more than a few meteors per 
hour. 



1.2 Dynamical Modeling of the Draconid Stream 

An important aspect of meteor physics is that it can re- 
veal physical properties of small solar system bodies with- 
out the cost of a spacecraft mission. Meteor streams are 
formed from ejections occurring during multiple perihelion 
passages of a parent body as well as occasionally through 
catastrophic break-ups. The meteoroids released during each 
ejection epoch will undergo a slightly different orbital evolu- 
tion compared to the parent body, depending on their mass 
and ejection state, but generally they are diluted into the 
background (annual) stream which is usually too sparse to 
allow a statistically sufficient sampling of shower meteoroids 
and typically includes meteoroids from many ejection eras. 
However, recently ejected material from the parent body 
may encounter the Earth and produce a meteor outburst, 
often comprising material from a single ejection epoch. Ob- 
servations of such outbursts allow us to directly sample a 
large number of meteoroids of known age from the parent 
body. The interaction of these meteoroids with the Earth's 
atmosphere can reveal clues about their properties, such as 
mass and chemical composition, which are directly related 
to the parent body. For example, the mismatch of visual and 
radar activity in the 2005 Draconid outbursts may suggest 
that the 2005 return was dominated by an abundance of 
sub-visual meteoroids from the 1946 perihelion pas sage of 
21P/Giacobini-Zinner (|Campbell-Brown et allbOOSl '). 

The dust trail theory of meteor outbursts/storms de- 
veloped in the 1990s turned out to be a robust way to un- 
derstand the linkage between outbursts of meteor showers 
and both the activity and ejection age from the p arent body 
(|Emerianen"kolll992l : lAshe3ll999l : lArlt et al.|[l99a '). However, 
to apply such a theory with confidence, sufficient observa- 
tions of the meteor shower itself and its parent body are 
needed as constraints. Virtually, no meteor showers can sat- 
isfy this requirement except the Leonids, mostly because 
of a lack of significant trail encounters and outbursts that 
might allow us to accumulate sufficient observations about 
individual dust trails. 

Prior to the late 1990's, potential Draconid outbursts 
were predicted by considering the distance between the or- 
bits of 21P/Giacobini-Zinner and the Earth, the time differ- 
ence between the arrival of the two bodies to the intersection 
point, and the geometry of comet's perihelion relative to the 
Earth's orbit. The predictions of Draconid outbursts in 1946, 
1985 and 1998 using this metho d were in good agreement 
with observations (Watson 1946i: iSpaldingj ll985: Langbroekl 
Il997h . while the predictions for 1952, 1959, 1972 and 1979 
did not agree. This reflects the value of applying the dust 
trail model to young, outburst streams such as the Dra- 
conids. 
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Figure 1. Location and geographic distribution of the main 
CMOR station (Zehr) and other remote sites as of October 2011. 



1.3 Outburst Predictions for 2011 



ISatd (|2003l ) first applied tiie dust trail theory to the Dra- 
conids and studied the trail encounter situ ations in 1998 and 
1999, while ICampbell-Brown et al. (200y) adopted a stream 
model proposed bv lVaubaillon et alJ (2005) to study the un- 
expected outburst of the Draconids in 2005. Further numer- 
ical simulations were carried out by a number of researchers 
(see Table [l] for a summary), and Earth's encounter of 
21P/Giacobini-Zinner's ejections from 1866 to 1907 (here- 
after the 1866-1907 "trails") were predicted to occur be- 
tween 16-20h UT, October 8, 2011. The associated peak 
visual rate p redictions ranged from 50 hr~^ ([Maslov 2011) 
to 600 hr"^ JVaubaillon et al.ll20lil ). Almost all observing 
techniques, including airborne missions, were deployed for 
the outburst, and preliminary analyses i ndicate that the 
outb urst arrived more or less as predicted (|Vaubaillon et al.l 

In this paper, we present the observation and analysis 
of the 2011 Draconid outburst as recorded by the Canadian 
Meteor Orbit Radar (CMOR). We focus on the following 
topics: 

(i) characteristics of the outburst, such as radiant, veloc- 
ity distribution, mass distribution, and flux variation; 

(ii) comparison with dynamical stream model predic- 
tions; 

(iii) meteoroid structure as revealed through ablation 
modeling. 



2 RADAR INSTRUMENTATION 

2.1 The CMOR System 

The Canadian Meteor Orbit Radar (CMOR) is an interfer- 
ometric backscatter meteor radar located near London, On- 
tario, Canada that is designed to observe meteor echoes and 
perform basic analysis continuously and automatically. It is 
based on the comme rcially available SKiYMET system (e.g. 
iHocking et al.ll200ll ) with some modificati ons to optimize for 
astronomical meteor echo detection (e.g. Ijones et al.ll2005l : 



Figure 2. Simplified example of how CMOR measures meteor 
trajectories. In this example, three radar sites detect signals re- 
flected from the meteor trail at different points as the meteoroid 
moves in the atmosphere. The time differences between the three 
observations, together with the interferometric direction mea- 
sured from the main site, can be used to construct the trajectory 
of the meteor. 



IWervk fc Brownll2012l ). Currently it consists of six sites (Fig- 
ure [U and operates at 29.85 MHz at 12 kW peak poweo 
The radar detects meteors through reflection of a transmit- 
ted pulse from the ionization trail left behind during meteor 
ablation and subsequently received after specular reflection 
by receivers. Meteors observed by the radar from the main 
station (Zehr) are always 90° from the apparent radiant. Be- 
cause of this geometry, the echo rate for radar observations 
has a secondary minimum when the radiant has a zenith 
angle <~ 30°, as opposed to visual, photographic and video 
meteor observation, where the minimum in apparent rates 
typically occurs when the radiant is near the horizon. 

If the echo from a meteor is detected at A*' sites (A'^ ^ 3), 
we can record A'' specular scattering positions along the me- 
teor trail, allowing us to measure the trajectory of the me- 
teor using_a__tinic-of-fliglit_(tof}_algorithni (Figure (2] see 
also I Jones et al . 2005; Wervk fc Br own 2012). However, the 
uncertainty in such a trajectory largely depends on echo 
strength and geometry. If the echo is weak, determination 
of the time of occurrence (or "time pick") will be difficult. 
Since the specular locations along a trail may coincide from 
different sites depending on the geometry of the meteor tra- 
jectory, the time difference between these sites may be small, 
and small uncertainties of the time pick may result in sig- 
nificant errors to the trajectory. 

Meteor echoes detected by CMOR are automatically 
processed to remove bad detections and to correlate com- 
mon eve nts across the sites . Details of th i s process ca n be 
found in I Jones et al.l (|2005| ): iBrown et all l|2008l . |2010| ) and 
IWervk fc Brownl (|2012l V 



^ The system also operates in 17.45 and 38.15 MHz, but only 
data at 29.85 MHz is used in this study. 
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Table 1. Summary of Draconid 2011 outburst predictions of time and radiant position for individual trail encounters. The forecast data 
from Lyytinen and Moser are quoted from http://draconids.seti.org/ (accessed Mar. 17, 2013). 



Trail 


Reference 


Predicted ag, 5g 
(J2000) 


Predicted maximum 
(Oct. 8, UT) 


1866 


Vaubaillon et al. (2011) 


263.3°, +55.3° 


16; 13 


1873 


Vaubaillon et al. (2011') 


263.2°, +55.4° 


16;29 


1880 


Vaubaillon et al. (2011) 
Sato r2003l 


263.2°, +55.4° 


16:53 
19:04 



1887 



Vaubaillon et al. (2011) 

.Sato (2003) 

Lyytinen 



263.2°, +55.4° 



1900 



Vaubaillon et al. (2011) 

Sato (2003) 

Moser 

Lyytinen 



263.2°, +55.8° 



17:25 
17:05 
17:02 





Maslov (2011) 


263.2°, 


+55.4° 


17:04 








1894 


Vaubaillon et al. (2011) 
Maslov (2011) 


263.2°, 
263.3°, 


+55.4° 
+55.5° 


18:45 
18:06 









20:01 
20:36 
19:52 
20:12 





Maslov (2011) 


263.3°, +55.8° 


20:42 


1907 


Vaubaillon et al. ("2011) 
.Sato (2003) 
Maslov (2011) 


262.5°, +55.4° 
263.2°, +55.8° 


19:26 
19:59 
21:05 




^^^^^^ 





Table 2. Basic specification of the 29.85 MHz CMOR system, adapted from lWervk fc Brownl l|2012l ). 



Parameter 



Value 



Frequency 


29.85 MHz 


Range interval 


15-255 km 


Range resolution 


3 km 


Pulse frequency 


532 Hz 


Peak power 


12 kW 


Noise floor 


-107 dBm 


Dynamic range 


33 dB 


Beam size 


55° at -3 dB point 



2.2 Echo Types 

CMOR's antennas have broad beams, allowing the radar to 
detect meteors appearing almost anywhere in the sky. The 
gain distribution resembles a bubble rather than a uniform 
distribution. This not only afTects the number of meteors de- 
tected with respect to different ranges, but also affects the 
determination of the physical properties of individual mete- 
oroids. One of the most significant properties of any radar 
echoes is the trail type of the meteor, namely underdense or 
overdense. 

An underdense trail occurs when radio waves scatter 
from all the individual the electrons in the trail. In contrast, 
an overdense echo occurs when the radio wave cannot com- 
pletely penetrate the meteor trail due to the trail plasma fre- 
quency being higher than the radar wave frequency. Strictly 
speaking, echoes with appearance between underdense and 
overdense should be classified as "transition" echoes, but for 
simplicity we also consider them as overdense echoes in this 
study. 

Visually, the amplitude-time series of an overdense echo 
will appeared as "fiat" (i.e. does not decay) for some time un- 



til ambipolar diffusion makes the trail underdense. Roughly 
speaking, an overdense echo is associated with a larger mete- 
oroid; but being overdense or not also depends on the specu- 
lar geometry to some degree. If an echo falls within a region 
where the antenna is weak, the chance of being overdense is 
higher. 

Examples of underdense and overdense echoes are 
shown in Figure |31 Detailed theory of these two echo types 
is beyond t he scope of this paper, but interested readers 
may refer to iMcKinlevI l|l96ll . §8) or ICeplecha et al.l (Il998l . 
§4) for details. The CMOR automatic detection algorithms 
are tuned to accept underdense and transition echoes but 
generally suppress overdense-type echoes. 



3 OBSERVATIONAL METHODOLOGY 

3.1 Selection of Draconid Echoes 

For our study, we first separated the Draconid meteor echoes 
from other echoes. To include as many Draconid meteors as 
possible, we use two methods to select underdense and over- 
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dense Draconid echoes respectively; this is briefly described 
below. 

For multi-station echoes, we developed user-interactive 
software to view, select, and analyze the geometry and tra- 
jectory of meteors automatically detected in the raw data 
stream. These are mainly underdense echoes with some 
transition-type echoes included as well. We consider any me- 
teor radiant within 5° and 20% of the predicted 1900-trail 
geocentric radiant, (i.e. Og = 263.2°, Sg = -1-55. 8°, J2000), 



and geocentric velocity Vg = 20.4 km-s~^ (IJenniskend 
120061 ). to be Draconid candidates. The velocity constraint 
is broad as Draconids are more fragile than other meteor 
streams IjJacchia et al.lll950l : lFiocco fc Colombdll964l '). lead- 
ing to larger uncertainty in measured velocity as a result of 
deceleration with height of detection. Using this method, we 
selected 61 Draconid candidates, including 39 underdense 
echoes and 22 transition-type echoes. We refer to these as 
the complete dataset. 

Because of the timing of the outburst (~ 20h UT) and 
the geometry of the radiant for CMOR at this time, we ex- 
pect many overdense echoes to be detectable from the Dra- 
conids. However, as some non-specular overdense echoes are 
detected through wind-twisting or development of plasma 
irregularities (Figure O), they can be difficult to detect and 
measure in an automated fashion. Therefore, we manually 
examined the raw data from the main site (Zehr) between 
Oct. 7-9, 2011 at 15-20h UT each day for overdense Dra- 
conid echoes. The trajectories of the trails were manually 
determined for cases when they were well-observed at other 
remote sites. This permitted us to isolate overdense Dra- 
conids. Trails with observations from less than three stations 
but with interferometric position consistent with Draconid 
meteor (i.e. at right angles to the Draconid radiant), were 
marked as "possible" Draconids. Following this method, we 
identified 148 overdense specular Draconid echoes as well as 
31 possible Draconid echoes, which formed another sample 
(Figure |4l referred to as the overdense dataset). For the rel- 
ative flux analysis, the number of possible overdense echoes 
on October 8, 2011 is subtracted from echoes from the same 
periods on Oct. 7 and 9, 2011, to correct for possible spo- 
radic overdense contamination. 

3.2 Uncertainties 

As the complete dataset was measured in an automatic fash- 
ion, it is possible to objectively estimate the uncertainty 
in the data. This involves an estimate of the uncertainty 
of the position of to (the specular point at each station) 
and the interferometric direction, both of which are directly 
measured, to estimate the error in the tof velocity. To do 
this, we used an echo simulator, which constructs a synthetic 
echo based on the observing geometry and strength of the 
mea sured signal assu ming a Gaussian distribution of noise 
fsee lWervk fc Brownll2012l '). Each pseudo-echo is then ana- 
lyzed using the same detection/measurement algorithm that 
is used for real observations. A comparison is then made be- 
tween the measured model results and the "true" (input) 
state to determine the expected random observational error 
in radiant and speed measurements, under the assumption 
of a perfect underdense specular echo. We ran the simula- 
tion 65 536 times for each echo and recorded the standard 
deviation to be the uncertainty in speed and radiant. We 
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Figure 3. A typical Draconid underdense echo (above), over- 
dense/transition echo (middle) and wind twisted overdense echo 
(below). We define these by the shape of their amplitude— time 



series (see §2.2). 



were able to use this technique on 41 echoes in the complete 
dataset; for the other 20 echoes in the complete dataset, we 
found it necessary to manually revise the time pick for re- 
liable results, so no uncertainty can be estimated for these 
echoes. We believe the error estimation for the 41 echoes 
is representative of the uncertainty of the entire complete 
dataset. 

The velocity of each echo in t he complete dataset is 
determined using the tof method (JBaggalev et al.Ml994l ). 
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Figure 4. Number-time series of definite and possible Draconid 
overdense echoes. 



To verify that the tof method yields accurate velocity 
measurements, we al so use the Fresnel phase-time method 
ICervera et al.1 119971 ) to measure the velocities of echoes 
in the complete dataset that exhibit such features. This 
method works as shown in Figure [S] when a specular meteor 
echo exhibits characteristic phase changes before the spec- 
ular (to) point, we can determine the speed of the meteor 
by fitting the slope of the phase change before the to point. 
Because the pre-to measurement precision does not depend 
on trail geometry with respect to radar sites (as opposed to 
the tof method), it can be an independent check for the lat- 
ter. We were able to manually measure the pre-to features of 
57 echoes in the complete dataset covering a wide range of 
heights (Figure [7|), and found they roughly agreed with the 
result from the tof -method where errors are largely due to 
geometry of tof and height averaging of the assumed con- 
stant speed (Figure [5|. We should note that the result from 
the pre-to method is only used to broadly verify the tof- 
method, as the pre-to Fresnel methoqj requires subjective 
determination of the starting point of phase changes, so the 
uncertainty estimations are also subjective. 
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Figure 5. An example of the pre-to phase-change exhibited by a 
Draconid meteor echo (upper). We remove the 2tt jumps and use 
the phase-distance relation (d oc ^/0) to determin e the to point 
and the velocity (lower) following the method of ICervera et al.l 
lll99'il ). In this example, the pre-to velocity is determined to be 
21.8 kms-i. 




20 22 24 26 28 

Tof-velocity (km-s^) 

Figure 6. Comparison between tof-derived velocities and pre- 
to-derived velocities. Dashed line is the 1:1 slope. 



^ There is another Fresnel meth od called the amplitude-time 
method (see lCeplecha et al.lll998l , §4.6.1); however only 7 echoes 
in our sample exhibited such a feature, and therefore we did not 
use this method. We note that the small number of Draconid 
echoes displaying Fresnel amplitude variations is another indica- 
tion of the widespread fragmentation behavior of the Draconids. 
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Figure 7. Height distribution of meteor echoes with pre-to deter- 
mined velocities. The average deceleration trend is clearly repre- 
sented by the change in speed with scattering height. The dashed 
vertic al line represents the expected out of atmosphere speed from 
Ijenniskena t2006iV 



Figure 8. The distribution of orbital elements a, e and i for the 
complete dataset. Shaded areas depict the range of respective 
orbital elements of the p seudo particles in mass r ange of 10^^" to 
10~^ kg as simulated by I Vaubaillon et al.l 1 120111 ) (assuming bulk 
density of 300 kg ■ m~'^). 



4 ANALYSES, RESULTS AND DISCUSSION 

4.1 Radiant and Orbits 

The weighted mean geocentric radiant for all echoes in the 
complete dataset was determined to be ag — 261.9° ±0.3°, 
Sg = +55.3° ± 0.3° (J2000), as shown in Tableland Fig- 
ure 1101 The errors shown here are standard Poisson error 
(1/viV). We also examined the meteors detected around 
the time of two predicted main peaks (see Table [l}, namely 
the expected peaks from the 1873-1894 and the 1900-1907 
trails. We examined subsets of the complete dataset that 
spanned these predicted trail passage times and found the 
weighted (inversely by error) geocentric radiants. Compar- 
ing these to the predictions summarized in Table [TJ the ra- 
diants determined from our observations match the model 
prediction reasonably well (~ 1°), and are comparable to 
the observati ons of Shigaraki middle and upper atmosphere 
(MU) radar (|Kero et al.ll2012l) (Figure [TOl we note that the 
MU results span a large period of time and thus represent a 
mixture of meteoroids from many trails) . 

The multi-station observations in the complete dataset 
also allow determination of the orbits. The distribution of 
a, e and i is shown as Figure |8] We also selected a set of 
echoes with Aa/a is below 0.2 (Aa was determined by the 
echo simulator as described in §3), defining a set of high 
quality multi-station orbits. The distribution of these orbits 
is shown in Figure O We also plot the range o f a, e and 
i of the simulated particles by IVaubaillon et al.l l|201ll ) as 
shaded bars in each graph. It can be seen that most echoes 
are located in a £ (2, 3) which is smaller than the a ~ 3.5 of 
the simulated particles, as well as the parent body. 

What is the cause of this difference in Vg and a be- 
tween observation and theoretical values? One possible ex- 
planation for the smaller observed a (and h ence smaller 
Vg) a s co mpared to modeling, as noted by iBrown et al.l 
l|2004 ) and ICampbell-Brown et al.1 l|2006l) , is an underesti- 
mation of the deceleration correction for Draconid echoes us- 
ing the standard deceleration correction applied to CMOR 
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Figure 9. The distribution of orbital elements a, e and i for the 
echoes in the complete dataset which Aa/a < 0.2. Shaded ar- 
eas depict the range of respective orbital elements of the pseudo 
particles in mass rang e of 10"^" to 10~^ kg as simulated by 
IVaubaillon et al.l 1I2OIII ) (assuming bulk density of 300 kg ■ m~^). 



echoes as a whole. Since the correction was derived from 
othe r showers which are no t as fragile as the Draconids 
(e.g. iFiocco fc Colombdll964l \ the out-of-atmosphere veloc- 
ity (foo) for Draconids may therefore be underestimated, re- 
sulting in a lower apparent initial velocity corresponding to 
a smaller a. As shown in Figure fTTI instead of the corrected 
11 00 showing no height trend, we observe a slowly increasing 
iioo with respect to height, in dicating a str o nger d eceleration 
to the meteoroids than the iBrown et al.l (|200J) algorithm 
predicts. 

However, from Table [3] we also notice that the observed 
Vg for Draconid echoes occurred during the times of the ex- 
pected arrival of the 1873-1894 and 1900-1907 trails are 
even lower than the overall Vg. A quick data check shows 
that the echoes detected during the passage of the two trails 
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Table 3. CMOR-observed weighted mean geocentric radiant and velocity for all detected Draconid meteors. Also shown are average 
radiants and speeds expected to be associated with the 1873-1894 and 1900-1907 trails based on the echo time of appearance in 
comparison to the model in the complete dataset. 



Time period 


Oig, Sg 

(J2000) 


■"9 

(km-s-1) 


N 


Annual radiant bv Jenniskens (2006) 
CMOR observation: Oct. 7-9 


264.1°, +57.6° 
261.9° ±0.3°, +55.3° ±0.3° 


20.4 
19.1 ±0.3 


61 


Vaubaillon's prediction for 1873-1894 trails 
CMOR observation: Oct. 8 16:19-18:55UT (1873-1894 trails) 


263.2°, +55.4° 
262.2° ±0.4°, +54.9° ±0.4° 


18.3 ±0.4 


32 


Vaubaillon's prediction for 1900-07 trails 
CMOR observation: Oct. 8 19:00-20:00UT (1900-1907 trails) 


262.8°, +55.6° 
262.5° ±0.6°, +55.3° ±0.8° 


17.0 ±0.5 
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Figure 10. CMOR-observed weighted mean geocentric radiant for the 1873-1894 trails (upper left), the 1900-1907 trails (upper right) 
and all echoes in the complete dataset (lower left). Uncertainty bars based on echosim results for individual radiants are plotted where 
applicable. Grey circles indicate a 5° circle within the predicted mean radiant of the model predicted respective trail or annual radiant. 
Sh aded grey dots in the up per two figures represent the radiant of pseudo particles in the mass range of 10"^" to 10~^ kg as simulated 
bv I Vaubaillon et al.l 1120111 ') (assuming a bulk density of 300 kg ■ m~^). The o verall comparison between observations of CMOR and MU 
radar JKero et al.l 2012h . as well as the prediction bv I Vaubaillon et al.l 1 1201 j) . is given in lower right of the figure. 



were roughly 10-15% slower than the those detected in the 
early hours of Oct 8 (the weighted mean of echoes detected 
before 12h UT is Vg = 19.6 ± 0.7 km • s"^). Although a dif- 
ference in physical properties between these two trails and 
the Draconid background is a possible explanation, we sug- 



gest that the contribution of a small radiant zenith angle (rj) 
(Figure I12p is more convincing. Due to a higher radiant el- 
evation, the meteoroid trajectory is more vertical, resulting 
in a deeper penetration through the atmosphere (Figure [TS] 



epf 



see also iVerniani 197S, Equation 25). We note that the Vg for 
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Figure 11. Height versus deceleration corrected velocity deter- 
mined using the tof method corrected for mean deceleration at 
infinity (tioo) for all detected echoes in the complete dataset. If 
the deceleration correction was correct, the data points should not 
show any trend with height. The vertical line represents the th e- 
oretical infinity velocity (23.3 km ■ s~^) from I JenniskensI 1120061) 



Figure 13. Specular height versus radiant zenith angle (r]) for 
all detected echoes in the complete dataset. A trend can be seen 
where the specular height is lower for smaller rj, supporting the 
argument that trails that were more vertical were able to pene- 
trate to lower heights. 
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Figure 12. Geocentric speed (vg) versus radiant zenith angle 
{rj) for all detected echoes in the complete dataset. A trend can 
be seen where Vg is lower for smaller rj, where trails were more 
vertical and were able to penetrate to lower heights. 



1900-1907 trails derived from European video observations 
ag rees the expected value {vg = 20.9 ± 1.0 km • s~^ as given 
bv lJenniskens et al.l201ll ) which supports our argument that 
the CMOR speeds are lower due to underestimation in de- 
celeration. 



4.2 Mass Distribution 

One of the diagnostic parameters one can get from statis- 
tical analysis of a meteor outburst is the mass distribution 
index s. The mass index is a measure of mass distribution 
of the meteoroids; it is defined such that the number of me- 
teoroids between mass interval m and m + dm is a power 
law described by m"". 



It is very difficult from specular radar echoes to directly 
measure m for each meteoroid; however, in radar theory, m 
and the electron line de nsity q are line arly related under 
certain assumptions fe.g. lMcKinlevlll96ll . §7.3-7.5). Consid- 
ering that q is proportional to echo amplitude for underdense 
echoes, we can simply substitute amplitude for m, and use 
the slope of cumulative logarithmic amplitude distrib ution 
to derive s (|McIntoshl [T96a : ISimek fc McIntoshlll968l '). For 
overde nse echoes, a di ffusion-limited echo duration implies 
T cc q (JMcKinlevI 1 196 ll . §8.9), and we can use echo du ration 
instead of amplitude to estimate s (|McIntoshl 1 19681 ) inde- 
pendently. 

We have few underdense echoes due to the poor radi- 
ant geometry at the outburst peak for CMOR, therefore we 
first use the echo duration method to estimate s. However, 
since overdense echoes can be affected by winds (Figure [3]), 
there is a risk t hat the overdense echoes may not be spec- 
ular (see [McKinlcy 1961, §4.9), leading to possible sporadic 
contamination. To address this issue, we manually inspected 
every echo for its pre-to phase feature (Figure^. Only spec- 
ular echoes exhibit such a feature, and therefore we can re- 
move non-specular echoes. In the end, we selected a total of 
155 Draconid overdense echoes on Oct. 8 in the overdense 
dataset. 

The identification of an overdense echo is rather sub- 
jective, leading to doubt as to whether some underdense 
echoes may have been confused with short-duration over- 
dense echoes. However, unlike the duration of overdense 
echoes, which depends on radar wavelength, ambipolar diffu- 
sion coefficient, and electron line density, the duration of un- 
derdense echoes only depends on the radar wavelength and 
a mbipolar diffusion coefficient (and hence height), as shown 
in lMcKinlevI l|l96ll . §8.5). Combining t his with the empirica l 
diffusion coefficient formula derived bv lMassev et al.l ([197l|), 
we have 

^ 16 ■ TT^ ■ 10006h-4.74 ^ ' 

where tud is the duration and A is the radar wavelength. 
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Figure 16. Mass indices for Draconid underdense echoes between 
Oct. 7-9 determined using the underdense amplitude method. 



Figure 14. The mass index of Draconids on Oct. 8 determined 
with selected echoes in the overdense dataset. The turnover from 
diffusion-limited to chemistry-limited duration can be seen at 
tc = 2.7 s. The mass distribution index determined by fitting 
the diffusion-limited portion is 1.72 ± 0.01. 
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Figure 15. The height range of the selected overdense echoes 
in overdense database. It can be seen that the lower limit of 
echo height is ~ 75 km. The shaded area marks the underdense 
region defined by Equation 1. We note that if this population 
were mainly underdense, a clear duration vs. height trend would 
be present. 



both in SI units, and h is the height in kilometers. With 
A ~ 10 m and a lower limit of /i ~ 75 km (Figure [TSJI. we see 
the upper limit of t^d to be ~ 1 s. Based on this result, it 
is reasonable to label any CMOR echo with duration longer 
than ~ 1 s as an overdense echo. 

We found s = 1.72±0.01 by fitting the linear portion of 
diffusion-limited overdense duration data poi nts with Fig - 
ure [Til and using the relation N ex r^fi"'')/'* (|Weisslll96]| ). 
where A'^ is the cumulative number of echoes. The uncer- 
tainty here only depicts the uncertainty of least-squares fit. 

We can also measure s using the underdense echo am- 
plitudes detected at the main site (Zehr). The advantage 
of this method is that it will include more events, includ- 
ing the possible Draconids that were filtered out in multi- 



station trajectory measurements. Sporadic contamination is 
possible, but given the low background rates at this time 
of the day (~ 15h local time), we expect it to be represen- 
tative of the mass index of the outburst. We select all the 
echoes perpendicular to within 5° of the apparent radiant 
and Vg = 20.4 km • s~^ (derived using rise-time estimation 
which can be applied to single station data) over a range 
of 10%. We find 84 single station echoes match our accep- 
tance criteria for the period of Oct. 7-9, and the mass index 
is determined by using the relation A^ ex A^'" (jMcIntoshl 
Il968 ) (where A*' is the cumulative number of echoes exceeds 
amplitude ^4). As shown in Figure [16] the mass index de- 
termined in this way for Oct. 8 was s — 1.75 ± 0.01. The 
uncertainty given here is from least-square fitting only and 
is several times smaller than the real unce rtainty given the 
small sample size (see iBlaauw et al.l 1201 ll . for discussion) . 
The mass indices for Oct. 7 and 9 appear to be higher than 
Oct. 8, possibly suggesting sporadic contamination, but they 
are subject to low confidence since the sample sizes are too 
small for quantitative estimates of s (14 for Oct. 7 and 16 
for Oct. 9). We are unable to compare this with the mass in- 
dices determined from overdense echoes for those two days, 
as there are also very few Draconid overdense echoes (< 10) 
for either of the days. 

The s — 1.7 5 on O c t. 8, the peak, is low compared 
to that found by ISimekl (|l994[ l (s = 2.06 for underdense 
echoes, s — 2.11 for overdense echoes) for the 1985 Draconid 
outburst. The implies that the 2011 outburst was richer in 
bright meteors than the 1985 outburst. Of course, such a 
conclusion is weakened by the small sample size that we 
used, but we note that it is consistent with the theoretical 
suggestion by Maslo^j, whose model suggested that the 1985 
outburst was composed of "quite faint meteors" while the 
2011 outburst was relatively rich in bright meteors. Com- 
bining the value derived from overdense and single station 
data, we suggest the 2011 outburst had s ~ 1.75, appropri- 
ate to -f 3 ^ Mv ^ +7. 



^ jhttp : //f eraj .narod.ru/Radlants/Predlctions/1901-2100eng/Draconidsl90~ 

(accessed Mar. 17, 2013). 
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4.3 Flux 

We counted the events in the overdense dataset and se- 
lected events around the time of the maximum in the effec- 
tive specular collecting area (17:30-18:30 UT) to compare to 
the zenith hourly rate (ZHR) published by the International 
Meteor Organizatioijj, shown in Figure [T71 As noted in pre- 
vious sections, the main peak measured with visual data for 
the 2011 outburst occurred under very poor specular scat- 
tering geometry for CMOR, such that only the outburst rise 
is measurable. However, the scaled ZHR from the overdense 
radar data shows a similar variation to the visual data. 



4.4 Meteoroid Structure Inferred from Ablation 
Modeling 

The meteoroid ablation mo del developed by 
ICampbell-Brown fc Koschnvl (|200J) was used to ex- 
plore the structure of radar-detected Draconid meteoroids. 
We can use the velocities and electron line densities mea- 
sured from different radar stations as constraints, and fit 
the model to match the observations. The velocity used 
here is the tof velocity described in §3; the electron line 
dens ity is computed from the amplitude/power of the echo 
(see ICeplecha et all Il998l . §4.2, §4.3). The computation 
of electron line density differs between underdense and 
overdense echoes. Following the definition of classic radar 
meteor theory, we use the underdense for mula for echoes 
with (J < 2.4 X 10" m"^ and vice versa (|McKinlevlll96ll . 
Equation 8-29). In this way, we perform the ablation entry 
modeling using two approaches: (i) modeling the Draconids 
as a "mean" population using all events; (ii) modeling 
individual events that showed noticeable deceleration across 
multiple stations. 



Table 4. Fixed 

ICampbell- Brown 



in put 

Koschnvl I2OO4 



parameters of the 

^__^^__^^^^^^^^ «a^^^, model. The infinity ve- 
locity Voo is from Ijenniskenj 1120061 ) ; the zenith angle uses the 
most probable value from o ur data; the bulk and g rain densities 
use the values suggested by iBorovicka et al.l 1120071 ) ; the heat of 
ablation and thermal conductivity u se the default values used by 
ICampbell-Brown fc Koschnvl | |2004|) . 



Parameter 



Value 



Infinity velocity, iioo 

Zenith angle, rj 

Bulk density, p^uifc 

Grain density, Pgrain 

Heat of ablation, q 

Thermal conductivity, k 



40° 

300 kg ■ m-3 

3 000 kg ■ m-3 

3 X 10*^ J ■ kg-i 

0.2 J-m-i -s-i -K- 



Table 5. Summary of parameters for the sensitivity test of 
fngrain and Ugrain for Figure [Tol and Figure l20l 



Figure 


irigrain (average) 


'pgrain 


Total mass 


Figure [T9l 


10"^ kg (fixed) 


20 

500 


2 X 10"^ kg 
5 X 10-5 kg 


Figure [20I 


2 X 10"* kg 
5 X lO-'^ kg 


100 (fixed) 


2 X 10-'^ kg 
5 X lO--"' kg 



density and velocity simulations, we see a trend that fixed- 
ngrain curves match the observations better; for the fixed- 
rrigrain curvcs, the sharp maximum of electron line densities 
near 93 km does not match the observed data points over 
the same height range. Also the velocity simulation showed 
little variation with different Ugrain in contrast to the obser- 
vation. 



4-4-^ Modeling Draconids as a Population 

In this approach, we mainly focused on two meteoroid 
structural parameters: grain number (ngrain) and grain 
mass {mgrain)', for other tunable parameters, we either 
use the known properties of Drac onids or those used by 
ICampbell-Brown fc Koschnvl (I2OOJ), as summarized in Ta- 
bled 

First, we find starting values for {mgrain, ngrain)- We fix 
the total meteoroid mass in the range [10~^ kg, 10~* kg] , and 
plot the results of three combinations with {mgrain, ngrain) 
ranging from (10"" kg, 1) to (10"^'' kg, 10 000) to compare 
with radar observations, as shown in Figure 1181 Our best 
fit is near mgrain £ [10"* kg, 10"® kg] and Ugrain = 100, 
though the solution is not unique and we emphasize that we 
are in effect trying to fit an average Draconid, when there 
may be significant variation within the stream. 

Next we want to see which of the two parameters 



(mg,. 



atn, "'grain 



dominates the process. We investigate this 



by fixing either mgrain (to lO""" kg) or Ugratn (to 100), let- 
ting the other parameter, as well as total meteoroid mass, 
being variable as summarized in Table (5] The results are 
shown in Figure [T^ and Figure !^ For both the electron line 



^ |http : //www ■ Imo ■ net/llve/draconlds201 1/ , (accessed Mar. 17, 
2013). 



4-4-^ Modeling Individual Draconid Echoes 

As the time and phase measurements in CMOR data are 
synchronized, it is possible to use the Fresnel amplitude 
fluctuation features just after the to point (Figure [ 2T|| to 
precisely measure the positions (see ICeplecha et al.l Il998l . 
§4.6.1 and references therein); we therefore searched in the 
complete dataset for any echo showing significant decelera- 
tion using this method. We have many echoes which show no 
measurable deceleration, but this does not provide a strong 
constraint on the entry model parameters, so we focus in- 
stead on events with clear deceleration. Since the same me- 
teor seen at different sites usually corresponds to different 
positions on the trail, if the Fresnel feature shows up in ob- 
servations from more than one site, we can obtain a series of 
contiguous position measurements along the trail. If there 
is no deceleration, then the position should increment lin- 
early with time; therefore, by plotting position-time series, 
we can readily see if there is any deceleration. This is anal- 
ogous to the dis t ance along trail versus time plots used by 
IBorovicka et al.l (12007J) for video data to reveal meteoroid 
deceleration. In this way, we identified one event (detected 
on 17:31:47 UT on Oct. 8 at the main site) that showed sig- 
nificant deceleration f Figure E^ . No other Draconid showed 
measurable deceleration. 

To model this event, we start from the same set of pa- 
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Figure 17. The flux variation of Draconid outburst on Oct. 8 as observed by CMOR. Complete stands for ttie number of echoes in 
the complete dataset, overdense - possible stands for the possible Draconid overdense echoes in the overdense dataset, and overdense - 
definite stands for the definite Draconid overdense echoes in the overdense dataset (see §3.1 for details). The overdense data is scaled 
relative to the visual data from IMO, and plotted on the same scale. 



rameters given in Table H] except the zenith angle and entry 
velocity are based on the specific echo observation (measured 
to be 36.8° and 24.3 km ■ s"'^ from the specular position and 
pre-to velocity measured at the main site). We take two ap- 
proaches: for the first, we fix the bulk and grain density to 
the one recommended by JBorovicka et al.l l|2007h (i.e. the 
one used in the previous section, ph^iih = 300 kg • m~'^ and 
pgrain = 3000 kg ■ m~^), then tune rugrain and Ugrain until a 
minimum difference between the modeling and the observa- 
tion is achieved. The second approach leaves both the bulk 
and the grain density as well as the zenith angle and entry 
velocity flexible to eliminate any remaining differences. The 
results found with the two modeling approaches, as well as 
the input parameters, are given in Figure [23] and Table [6] 



4-4-3 Comparison and Discussion 

We compare o ur observational and modeling result to 
JBorovicka et al.l ([2007), who reported six decelerated Dra- 
conid meteors with sizes comparable to radar meteors, de- 
tected by double-station video system during the 2005 Dra- 
conid outburst. Their grain masses were in the range of 
10"" to 10"^ kg, with ptuik = 300 kg • m"^ with the as- 
sumption of Pgrair ~ "^ """ ^'" """^ 

grains between 10 
significant deceleration with turnover height (i.e. the height 
that V becomes significant) around 95 km. The differences 
with our results include a much smaller grain number (10^) 



ratn ~ 3 000 kg • m , and total number of 
'** and 10®. These six meteors all showed 



and a much larger grain mass (at the order of 10~^ kg) de- 
spite the density configurations being similar. Our 17:31:47 
event also showed a slightly lower turnover height (~ 93 km) 
than the video meteors. 



How to explain these differences? One possibility is 
the outbursts from 2005 and 2011 were ejected by Comet 
Giacobini-Zinner in different years: the material encoun- 
tered in 2011 was ejected earlier than 2005 (1873-1907 ver- 
sus 1 946, per lCampbell-Brown et aklbood : IVaubaillon et al.l 
I2OIH ). However, we note that it is still difficult at this stage 
to uniquely define the physical properties of meteoroids in 
this way due to the lack of observational constraints. For 
example, the statistical model fits shown in Figure [18] are 
broad averages which allow much uncertainties; for the mod- 
eling of the 17:31:47 event we have shown that a reasonable 
fit can be found with different sets of parameters, and sim- 
ilar phenomenon have been noted by the Borovicka et al. 
(see §7.1 of their paper). All these factors cast doubts on 
whether we are indeed seeing the physical differences be- 
tween meteoroids or just fit variations admittedable within 
the uncertainties. However, we still see that it is possible to 
use the observations to broadly constrain meteoroid proper- 
ties to some extent, and we feel encouraged by the fact that 
o ur results have shown b asic agreement to the earlier work 
of JBorovicka et all (|2007h . 
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Figure 18. Modeling results of three {mgraim "n-grain) coinbinations together with radar observations. The total metcoroid mass is fixed 
at the range of [10~® kg, lO^** kg], while the combination of {nigraini'n-grain) varies from ([10~® kg, 10"'* kg], 1) for the upper set, 
([10~ kg, 10~ kg], 100) for the middle set, and ([10^ kg, 10~ kg], 10 000) for the lower set. The curves stand for the modeling result 
using certain {nigrain, ingrain)- It can be seen that the best fit is the middle set, with nigrain G [10^* kg, 10^® kg] and rigrain = 100. 
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Table 6. Input parameters of two approaches used to model the ablation of the 17:31:47 event. For the 1"* approach, the bulk and grain 
densities arc set to tose given by|Borovicka ct al. (2007), while for the 2""^ approach the two are fit by the model. 



Parameter 



I''* approach 



approach 



Number of grains, rigrain 

Grain mass (average), rrigrain (kg) 

Grain mass range (kg) 

Total mass, m (kg) 

Infinity velocity, V(x, (km ■ s~^) 

Zenith angle, ri 

Bulk density, pf,„;fe (kg ■ m"^) 

Grain density, pgrain (kg • m"^) 
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Figure 19. Modeling results with fixed rrig 



10 kg and 



variable n„ 



20, 500, corresponding to total masses of 



2 X IQ- kg and 5 X lO-"" kg. Size and color of data points in 
velocity figure indicates number of sites: darker and larger points 
indicates more sites and vice versa. The vertical line for each data 
point depicts the height range constrained by multi-site observa- 
tions. Curves depict modeling solutions under different rigrain- 
The modelled fit does not show as much scatter as the observa- 
tions compared with Figure [20l 



Figure 20. Modeling results with fixed rtg 



100 and variable 



iTigrain = 5 X IQ- ^ kg, 2 X IQ- ** kg. Corresponding to total masses 
of 2 X 10—® kg and 5 X IQ- ^ kg. Size and color of data points in 
velocity figure indicates number of sites: darker and larger points 
indicates more sites and vice versa. The vertical line for each data 
point depicts the height range constrained by multi-site observa- 
tions. Curves depict modeling solutions under different rigrain- 
Comparing with Figure [T9l we see a better match to the obser- 
vations especially at middle and lower heights. 



5 CONCLUSIONS 

In this study, we have analyzed CMOR observations of the 
2011 Draconid outburst, including 61 multi-station specular 
echoes and 179 single-station overdense echoes. Our main 
results are: 

(i) The radiant of the outburst was determined to be Qg = 
261.9° ± 0.3°, Sg = -f55.3° ± 0.3° (J2000) with the multi- 



station echoes, which agreed with earlier modeling forecasts 
to ~ 1°. 

(ii) The averaged velocity determined from the CMOR 
data was by ~ 10 — 15% smaller than the expected value 
(17.0 - 19.1 km ■ s-i versus 20.4 km ■ s-^), likely due to the 
undercorrection of deceleration, combined with the effect of 
radiant geometry during the peak hours; 

(iii) The mass distribution index determined from 155 
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Figure 21. Example of an underdense Draconid echo witli four 
Presnel maxima amplitude features noted (marked by arrows). 
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Figure 22. The deceleration of the 17:31:47 event, plotted in 
position offset versus height. We linearly extrapolate the linear 
portion of the data points (from a height of 96 to 92.5 km), the 
position offset is the distance between the observed position and 
the linear extrapolation. 



overdense echoes was s = 1.72 ± 0.01 assuming the echoes 
used for such determination were diffusion-dominated. Al- 
ternatively, we selected 54 possible Draconids in the sample 
of single-station underdense echoes and determined the mass 
index to be around 1.75. Combining these values, we suggest 
that the mass index of the 2011 outburst to be s =1.75 for 
the magnitude range of -1-3 ^ AIv ^ +7; 

(iv) We compared the counts of overdense echoes to the 
visual data, and a clear consistency can be noted. Unfor- 
tunately, the peak time of the outburst occurred with very 
poor scattering geometry for CMOR, such that no compar- 
ison can be made at the peak time. However the rise of 
the meteor rates between the radar and visual results were 
highly comparable; 

(v) We used the meteoroid ablat ion model developed by 
ICampbell-Brown fc Koschnvl (|2004l ) to explore the structure 
of decent radar Draconid meteoroids. Assuming bulk and 
grain densities to be 300 kg • m~^ and 3 000 kg • m~^ re- 
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Figure 23. The modeling result for the 17:31:47 event: upper 
plot — modeling and observed electron line density; lower plot — 
modeling and observed theoretical minus observational position 
assuming a fixed speed. The theoretical fit was extrapolated using 
the position from 96 to 92.5 km as reference. 



spectively, the model seems to infer that the grain number 
of CMOR-observed Draconid meteoroids is ~ 100 regardless 
of meteoroid sizes; 

(vi) We also identified a Draconid meteor that showed 
clear deceleration. Two modeling approaches were at- 
tempted, suggested the grain numbers to be ~ 1 000, grain 
mass between 10"^'^ to 10~^ kg, and the total mass to be 
around 2 x 10~® kg, which were in rough agreeme nt with 
the video results reported bv lBorovicka et al.l (120071 ). 

As this paper was being finished, an unexpected and 
much more intense Draconid outburst was detected by 
CMOR on Oct. 8, 2012. Very preliminary analysis suggests a 
spectacular ZHR at the order of 8 000, far above the storm 
level. In-depth investigation of the 2012 event will be ad- 
dressed in a separate paper. 
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